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Abstract 



We study the problem of estimating parameters of the limiting equation 
of a multiscale diffusion in the case of averaging and homogenization, given 
data from the corresponding multiscale system. First, we review some recent 
results that make use of the maximum likelihood of the limiting equation. In 
particular, it has been shown that in the averaging case, the MLE will be 
asymptotically consistent in the limit while in the homogenization case, the 
MLE will be asymptotically consistent only if we subsample the data. Then, 
we focus on the problem of estimating the diffusion coefficient. We suggest a 
novel approach that makes use of the total p-variation, as defined in [11] and 
■»^J- \ avoids the subsampling step. The method is applied to a multiscale OU process. 
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1 Introduction 

It is often the case that the most accurate models for physical systems are large in 
dimension and multiscale in nature. One of the main tasks for applied mathematicians 
is to find coarse-grained models of smaller dimension that can effectively describe the 
dynamics of the system and are efficient to use (see, for example [121 EES 13 IE])- Once 
such a model is chosen, its free parameters are estimated by fitting the model to the 
existing data. Here, we study the challenges of this statistical estimation problem, in 
particular for the case where the coarse-grained model is a diffusion. Apart from the 
usual challenges of parameter estimation for diffusions, an additional problem that 
needs to be addressed in this setting is that of the mismatch between the full multiscale 
model that generated the data and the coarse-grained model that is fitted to the data. 
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A first discussion of this issue, in the context of averaging and homogenization for 
multiscale diffusions, can be found in [T5] EH EE] • 

A similar statistical estimation problem arises in the context of "equation-free" 
modeling. In this case, coarse-grained equations exist only locally and are locally 
fitted to the data. The main idea of "equation-free" modeling is to use these locally 
fitted coarse-grained equations in combination with a global algorithm (for example, 
Newton-Raphson) in order to answer questions about the global dynamics of the 
coarse-grained model (for example, finding the roots of the drift). In this process, 
we go through the following steps: we simulate short paths of the system for given 
initial conditions. These are used to locally estimate the effective dynamics. Then, we 
carefully choose the initial conditions for the following simulations so that we reach an 
answer to whatever question we set on the global dynamics of the system, as quickly 
and efficiently as possible (see [TU])- The statistical inference problem is similar to 
the one before: we have the data coming from the full model, we have a model for the 
effective local dynamics and we want to fit the data to this model. However, there 
is also an important difference: the available data is short paths of the full model. 
This issue has not been addressed in [15], [TJ] or P, where it is assumed that the 
time horizon is either fixed or goes to infinity at a certain rate. We will address this 
problem in section [3j by letting the time horizon T be of order O (e a ), where e is the 
scale separation variable and a > 0. Another important issue that we will address 
here is that of estimating the scale separation variable e. 

We will focus on a very simple Ornstein-Uhlenbeck model whose effective dynamics 
can be described by a scaled Brownian motion. This will allow us to perform precise 
computations, reach definite conclusions and build our intuition about the behavior 
of more general diffusions. We will only tackle the homogenization case and our goal 
will be to estimate the diffusion coefficient of the effective dynamics. This problem 
has also been addressed in [15], [1] . In both these papers, the diffusion coefficient is 
constant. In fact, in [1] the authors also focus on the Ornstein-Uhlenbeck model. Our 
main contribution is to demonstrate that in order to compute the diffusion coefficient, 
one should not use the quadratic variation commonly defined as a limit where we let 
the size of a partition go to zero but rather as a supremum over all partitions. This 
definition is discussed in [TT] and is at the core of the theory of rough paths, as it 
gives rise to a topology with respect to which the Ito map is continuous. 

In section [2] we review some of the core results for multiscale diffusions and their 
coarse-grained models. Then, we will review the results of [T5l [T^] and [I]. Finally, 
We will give a more precise description of "equation-free" modeling. 

In section [31 we go on to define a new set of estimators for the diffusion parameter 
of the coarse-grained model, in the case of homogenization. We perform explicit 
computations of their L2-error, which allows us to attest their performance. We 
conclude that they outperform the subsampled quadratic variance estimate studied 
in [15] Q]. Finally, we describe a heuristic way of estimating the scale separation 
parameter e. 
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2 MLE for multiscale diffusions: A review 

In this section, we review some of the main concepts that come into play in multiscale 
modeling. First, we describe the limiting equations for multiscale stochastic differen- 
tial equations. These allow us to reduce the dimension of the model. Then, we discuss 
the problem of the statistical estimation of parameters of the limiting equation given 
multiscale data and how this mismatch between model and data affects the result. 
Finally, we discuss a numerical algorithm that is applied when the limiting equations 
are completely unknown, which comes under the name of "equation-free" modeling. 



2.1 Limiting equations for multiscale diffusions 

Multiscale diffusions are a combination of two basic types of multiscale stochastic 
differential equations. The first is described by the following equations 

dX t = h(X t ,Y t )dt + ^{XuY^dWt 

(1) 

dY t = ±f 2 (X t ,Y t )dt+±a 2 (X t ,Y t )dVt 

where X t G X and Y t G y and X,y are finite dimensional Banach spaces. We call 
X the slow variable, Y the fast variable and e the scale separation parameter. The 
main assumptions are the following: 

Assumption 1. (i) The solution of the system exists, 
(ii) The equation 

dY x = \f 2 (x, Y t x )dt + -a 2 (x, Y t x )dV t 
is ergodic with unique invariant measure \x x , for every x G X . 
We expect that by the time X takes a small step A ~ O (1), 

f-t+A 

A 



^j t+ fi(X s ,Y s )ds* J h{XuVW t {dv) 



as a result of the ergodicity of Y . Similarly, 

^ J ai{X s ,Y s )ai{X s ,Y s )>ds w J a 1 (X t ,y)a 1 (X t , y y f i Xt (dy). 
where by (•)' we denote the transpose of a vector. We set 



AO) = J fi(?,y)(i x (dy), <ti(x)= a 1 (x,y)a 1 (x,y)'iJ lx {dy) 
and 

dX t = h{X t )dt + a^X^dWt. (2) 

We call (T5]) the averaged limiting equation and we call X the averaged limit. We 
expect that X t ~ X t , provided that they have the same initial conditions. Indeed, 
the following holds 
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Theorem 2.1 ([14]). Let X = T e and y = T d ~ £ . We assume that all coefficients 
in (QP are smooth in both x and y and that the matrix S 2 (x,y) = o~2{x,y)a2(x,y)' is 
positive definite, uniformly in x and y. Also, there exists a constant C > such that 

{z,B(x,y)z} >C\z\ 2 , \/(x,y)eXxy and z6l" 

where (•,•) denotes the Euclidean inner product. Then, if Xq = X , 

1^1 in C([0,T],X). 

Different types of convergence have also been proven under different assumptions 
(see [16] and [5]). 

The second basic type of multiscale stochastic differential equation is described by 
the following equations 

dX t = \h{X u Y t )dt 

(3) 

dY t = ±f 2 (X t ,Y t )dt+]a 2 (X t ,Y t )dV t 

where X t e X and Y t G y and X, y are finite dimensional Banach spaces. As before, 
we call X the slow variable, Y the fast variable. In addition to assumption [TJ we 
assume that 

Assumption 2. 

/ fi(x,y)/J>x(dy) = 0, Wx e X 
Jy 

where fi x as defined in assumption^ 

Then, we expect that by the time X takes a small step A ~ O (1), 
I rt+A i rt+A 



h{Xs,Y a )d 8 « — J h{X tl Y^)ds 

It follows from the Central Limit Theorem for ergodic Markov Processes (see [3]) that 
this will converge to a random number. More precisely, let us set 

p /*oo 

J x Jo 

and 

f(x) = (^2 J J™ f x {x,y) {P s h{x,-)){y)'li x {dy) 
where Pt are the transition kernels of the diffusion Y x . Finally, we set 

dx t = MX^dt + r(X t )dw t . (4) 

We call X the homogenized limiting equation. As before, we expect that X t ~ X t , 
provided that they have the same initial conditions. Indeed, similar to the averaging 
case, we can prove the following: 
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Theorem 2.2 ([H]). Let X = T e and y = T d ~ £ . We assume that all coefficients 
in (T3|) are smooth in both x and y and that the matrix S 2 (x,y) = o~2{x,y)a2(x,y)' is 
positive definite, uniformly in x and y. Also, there exists a constant C > such that 

(z,B(x,y)z) > C\z\ 2 , \/(x,y)eXxy ancb 6 R d ~ e , 

where (-,-) denotes the Euclidean inner product. Then, if Xq = X and assumption 
[H holds, we get that 

1^1 in C([0,T),X). 

Again, different types of convergence have also been proven under different assump- 
tions (see [16] and [5]). 

Theorems 12.11 and 12.21 allow us to replace the (X t , Y t ) system by X t . If we are only 
interested in the slow dynamics of the system, this allows us to reduce the dimension 
of the problem. For example, using the limiting equations we can simulate the slow 
dynamics of the process much faster, not only because of the dimension reduction but 
also because the dynamics of X do not depend on e. Thus, the step of any numerical 
algorithm used to simulate the dynamics can be of order O (1) rather than O (e 2 ) 
which would have been the case if we wanted to simulate the full multiscale system. 

2.2 Parameter estimation for multiscale diffusions: a review 

The theory reviewed in section 12.11 allows us to reduce the dimension of a multiscale 
system, approximating the slow dynamics by an diffusion of smaller dimension that 
does not have a multiscale structure anymore. In addition to multiscale diffusions, 
similar results hold for ordinary and partial differential equations (see [IS]). 

It is often the case that the dynamics of the full multiscale system - and conse- 
quently those of the limiting system - are not completely known. For example, in 
the case of multiscale diffusions, the drift and variance of the full system and thus 
the limiting system might depend on unknown parameters. This poses a statistical 
problem: how can we estimate these parameters give the multiscale data? In fact, it 
is even more realistic to ask to find the drift and diffusion coefficient of X given only 
X. This problem has been discussed in [151 EJ EE]- 

More precisely, in [15], the authors discuss the case where the drift of the limiting 
equation depends linearly on the unknown parameter while the diffusion parameter is 
constant. In [2], the authors extended the results of for generic drift but did not 
discuss the problem of estimating the diffusion parameter. Finally, in [T], the authors 
extend the results in [15] by also proving the asymptotic normality of the estimators, 
but the limit their study to the Ornstein-Uhlenbeck system. The approach taken so 
far is the following: 

(i) We pretend that the data comes from the limiting equation and we write down 
the corresponding maximum likelihood estimate (MLE) for the unknown pa- 
rameters; 

(ii) we study whether the mismatch between model and data leads to errors and, if 
so, we try to find a way to correct them. It has been shown that in the limit as 
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the scale separation parameter e — > 0, the MLE corresponding to the averaged 
equation is consistent. However, this is not true in the case of homogenization. 
The method used so far to correct this problem has been that of subsampling 
the data by a parameter 5. Then, for 5 ~ O (e a ) and a G [0, 2], it has been 
shown that the MLE that corresponds to the homogenized equation will be 
consistent in the limit e — > 0. Also, an effort has been made to identify the 
optimal subsampling rate, i.e. the optimal a. However, since e is usually an 
unknown, this is of little practical value. 

Note that a separate issue is that of writing the maximum likelihood of the limiting 
diffusion, which in the general multi-dimensional case can still be challenging (see 
j2j[9]). We will not discuss this issue here, however. 

We summarize the main results for the parameter estimation of the limiting equa- 
tions of multiscale diffusions in the following theorems: 

Theorem 2.3 (Drift estimation, averaging problem). Suppose that fi in |1)) depends 
on unknown parameters 9, i.e. f\(x) = fi(x;9). Let 9(x;T) be the MLE of 9 cor- 
responding to equation (Hj). Suppose that we observe {X t ,t G [0,T]} of system (Q]j 
corresponding to 8 = 9q. Then, under appropriate assumptions described in [L$ 
(Theorem 3.11), it is possible to show that 



lim dist (j)(X;T),9^j = 0, in probability 



where dist (•, •) is the asymmetric Hausdorff semi-distance and 9 t is a subset of the 
parameter space identified in the proof. Also 

limd# (9 e , 9 ) = 0, in probability 

e->0 

where dn (•, •) is the Hausdorff distance. 

Theorem 2.4 (Drift estimation, homogenization problem). Suppose that fx in ^ 
depends on unknown parameters 9, i.e. f\(x) = fi(x;9). Let 9(x;N,5) be the max- 
imizer of the discretized likelihood corresponding to equation (TJj] with step 5, where 
T = N5. Suppose that we observe {X t ,t G [0,T]} of system (TJ)) corresponding to 
9 = 9q. Then, under appropriate assumptions described in [Lffl (Theorem 4-5) and 
for 5 = e a with a G (0, 2) and N = [e -7 ] for 7 > a, it is possible to show that 

lim6'(X; N, 5) = 0, in probability. 

e-S>0 

The next two theorems deal with the estimation of the diffusion parameter of the 
limiting equation, given that this is constant. In that case, the MLE is the Quadratic 
Variation of the process. They assume that the dimension of the slow variable is 1. 

Theorem 2.5 (Diffusion estimation, averaging problem). Let X be the solution of (Q]) 
for a\ = 9 a constant. Then, under appropriate conditions described in fl5f (Theorem 
3.4) and for every e> 0, we have that 

m S ~ x ns\ 2 = 9 2 a.s. 



5^0 N5 

n=0 



where T = N5 is fixed. 
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Theorem 2.6 (Diffusion estimation, averaging problem). Let X be the solution of $3j) 
for f = 9 a constant. Then, under appropriate conditions described in f73] / (Theorem 
3.5) and for 5 = e a with a G (0, 1), we have that 

1 N-l 

&o m ^ l X ( n+1 ) 5 ~ X ns\ 2 = 2 a.s. 

n=0 

where T = NS is fixed. 

It is conjecture that Theorem 12.61 should hold for any a G (0, 2) and that the 
optimal a, i.e. the one that minimizes the error, is a — |. 

Clearly, the most interesting case is that of estimating the diffusion parameter of 
the homogenized system. This is the case that we will study in detail in section [31 
assuming that the process is an Ornstein-Uhlenbeck process. Also, note that when 
estimating the diffusion parameter, the length of the time interval T is fixed. We will 
relax this condition later on, for reasons explained in the following section. 

2.3 Equation-free modeling 

In practical applications it is often the case that the limiting equations fl2]) and fl4]) 
are completely unknown. More generally, let us say that we have good reasons to 
believe that a certain variable of a multiscale system that evolves slowly behaves 
like a diffusion at a certain scale but we have complete ignorance of its drift and 
diffusion coefficients. We would like to find a way to estimate these coefficients. In 
statistical terms, let us say that we are interested in the non-parametric estimation 
of the drift and diffusion coefficients of the limiting equation. Note that our data 
comes "on demand" but for a certain cost, by simulating the multiscale model for 
given conditions. 

A general algorithm for answering questions regarding the limiting dynamics of 
a quantity coming from a multiscale system that evolves slowly, when these are not 
explicitly known, comes under the name of "equation-free" algorithm (see [10j). In our 
case, this would suggest pairing the problem of local estimation with an interpolation 
algorithm in order to estimate the drift and diffusion functions, denoted by f(x) 
and &{x) respectively. We make this more concrete by describing the corresponding 
algorithm: 

0. Choose some initial condition xq and approximate f(x) and cr(x) by a local 
(polynomial) approximation around Xq. Simulate short paths of the multiscale 
system, so that the local approximation is acceptable. Note that the smaller 
the path, the better or simpler the local approximation. 

1. For n > 1, choose another starting point x n using the knowledge of f(x n -\) and 
a(x n -\) and possibly some of their derivatives on x n -i, according to the rules 
of your interpolation algorithm. 

2. Repeat step 0, replacing xo by x n . 

As mentioned above, the size of the path T needs to be small and possibly comparable 
to e. This is what led us to consider the estimation problem for T = e a . 
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3 The p-variation estimate 

In this section, we study the problem of estimating the diffusion parameter of the 
homogenization limit of a simple multiscale Ornstein-Uhlenbeck process. We hope 
that the detailed analysis will provide some intuition for the general problem. 
Consider the following system: 

dY^ = ^Y t 2 ' e dt 

(5) 

dY t 2 ' € = -j,Y^dt+\dW t 

with initial conditions Y ' e = y\ and Yq ,€ = y 2 . It is not hard to see that the 
homogenization limit as e — > is 

and the convergence holds in a strong sense: 

sup \Y^ e - yi -aW t \ ^0, ase^O. (6) 

*e[o,T] 

Note that for this particular example, Y t 1,e is exactly equal to 

Y?> e = Vl + aW t - ea (>f £ - y 2 ) 
and thus proving §6§ is equivalent to proving that 

e sup \Y 2 ' e — y 2 \ — ^ 0, as e — > 0. 

te[o,T] 

This follows from [6]. 

We want to estimate the diffusion parameter a given a path {Y^' e {u) ; t £ [0, T]}. 
If we were to follow the approach discussed in the previous section, we would use 
the maximum likelihood estimate that corresponds to the limiting equation. In this 
case, this would be the quadratic variation. However, as discussed earlier, this is not 
a good estimate since the quadratic variation for any fixed e > is zero. To correct 
this, we subsample the data, which leads to the following estimate: 

N 



NS 



i=i 

The asymptotic behavior of this estimate has been studied in [151 !]• I n i&ct, taking 
advantage of the simplicity of the model, we can compute the L 2 -error exactly, as a 
function of 6, e and N. We find that 




E 4 



e 4 _ * \ 2 3 + e"^ \ ( 1 
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For reasons explained earlier, we are interested in the behavior of this error not only 
when T is fixed but also for T — > 0. Thus, we set T = e a and, as before, 5 = e Q+/3 , 
which lead to N = e _/3 . We are interested in the behavior of the error as e — > 0. For 
these choices of T and 5, the square error will be 

E {a 2 -a 2 ) 2 ~0 ( e 4 ~ 2 (<*+« + e 2 ~ a + e?) (9) 

For a fixed, we see that the error will be small if < j3 < 2 — a. In fact, the optimal 
choice for j3 is /3 = i n which case the error becomes 

(E(a 2 -a 2 )y ~0( e V) (10) 

So, for a = 0, we get that the optimal sub sampling rate is f3 = |, which results to an 
optimal error of order O {^^j ■ However, if a > 0, the error can increase significantly, 

especially for non-optimal choices of 5. 

In the rest of this section, we are going to investigate the behavior of the p- variation 
norm as an estimator of a. The intuition comes from the following observation: we 
know that at scale 0(1), {Y t ' £ (u)) ;t G [0, T]} behaves like scaled Brownian motion 
while at scale O (e), it is a process of bounded variation (finite length). Could it be 
that at scale O (e a ), the process behaves like a process of finite p- variation, for some 
p that depends on a? If so, would the p- variation norm be a better estimator of cr? 



3.1 The total ^-variation 

We say that a real- valued continuous path X : [0, T] — > R has finite total p- variation 
if 

D p (X) T := sup [ V \X t(+1 -X tt \A <+oo, (11) 
\t e e-D([o,T]) J 

where T> ([0, T]) goes through the set of all finite partitions of the interval [0, T] (see 
also PH]). It is clear by the definition that a process of bounded variation will always 
have finite total p- variation for any p > 1. Also, note that the total p- variation as 
defined above will only be zero if the process is constant. Thus, the total p- variation 
of a non-constant bounded variation process will always be a positive number. 

For e > fixed, the process Y l,€ : [0, T] — > R defined in ([5]) is clearly of bounded 
variation, but its total variation is of order O (-). We will say that at scale O (e a ), 
the process Y 1,e behaves like a process of finite total p-variation in the limit if 

lim (D p (Y^) ea ) < +oo and Vg < p, lim [D q (Y x ' e )J = +oo. (12) 
We will prove the following: 

Theorem 3.1. At scale O (e a ) and 1 < a < 2, the process Y 1 ' 6 : [0,T] R defined 
in (TJJ) behaves like a process of finite total (2 — a) -variation in the limit. 



10 



First, we prove the following lemma: 

Lemma 3.2. Let X : [0,T] — > K be a real-valued differentiable path of bounded 
variation. Then, its total p-variation is given by 



D p (X) T := sup J2 ~ X t e \ P ) , (13) 

^[°' r D \t e ee([o,T]) J 

where £ ([0,T]) goes through all finite sets of extremals of X in the interval [0,T]. 

Proof. Consider the function 

f a , b (t) = \X t -X a \ p + \X b -X t \ p , a<t<b. 

This is maximized for t an extremal point (X t = 0) or at t = a or t = b. Thus, 
if D = {0, ti, . . . , t n -i> t n = T}, there exists a set of extremals £ with cardinality 
\£\ < n + 1, such that 

\ x t e +i - x t t \ p < \ Xt w - x tt\ p - 

t e ev t e e£ 

The set £ can be constructed by choosing T\ so that fo,t 2 {t) is maximized and so 
that f Tk _ 1 ,t k+1 {t) is maximized, for k — 2, . . . , n — 1. Thus 

sup f \ x t t+1 ~ XtA < sup ( \ x t e+1 -Xtf 

V ^ \t t €V([0,T\) J £ W> T » \t e e£(lo,T\) 

The opposite inequality is obvious and completes the proof. □ 

To prove the theorem, first we notice that 

D p (Y^) T = eaD p (Z 1 )r, (14) 

where (Z 1 ,^ 2 ) satisfy 

dZ] = Z 2 t dt 

dZ 2 = -Z 2 dt + dW t 

Now, Z l is clearly differentiable and thus, by the lemma 



d p (z%= sup Yl K +1 -K\ p 



£([0,T]) 



.t e ee([o,T]) 



The derivative of Z 1 is equal to Z 2 , so all its extremal points correspond to zero- 
crossings of Z 2 . So, for s,t e £, 

Z\ - Z\ = (W t - W s ) - (Z 2 - z 2 ) =w t -w. 
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and D p (Z 1 )t becomes 

D p (Z l ) T = sup V \ Wtl+1 -W te 
£ W> T » \t e es { [o,T]) 

= SI E i^w-^M , (15) 

~" \t e ee s ([o,T]) J 

where 

^([0,T]) = {0 = t ,ti,...,^(T),T} 

and {ti, . . . , tN s (T)} is the set of all zero-crossings of Z 2 in [0, T] that are at least 
distance <5 apart from each other, i.e. if tk G £5 ([O,? 1 ]) and & < N$(T), then is 
the first time that Z 2 crosses zero after time tk + 6. Note that the set of zero-crossings 
of Z 2 in [0, T] is an uncountable set that contains no intervals with probability 1. 
Equation (|T5|) follows from the following two facts: (i) in general, adding any point 
to the partition will increase the L p norm and thus the supremum is achieved for a 
countable set of zero-crossings and (ii) any countable set that is dense in the set of 
all zero-crossings will give the same result. 

If t§ is the stopping time of the first zero crossing of Z 2 after 5 given Zq = 0, then 
the random variables {r^ = (tk — t k _i),tk G £s([0,T]) ,k < N$(T)} are i.i.d. with 
the same law as that of t$. Thus, the sum J2t e ee s ([o t]) \ Wt i+1 — W te \ p is a sum of i.i.d. 
random variables of finite mean (to be computed in the following section) and as a 
consequence of the Law of Large Numbers, it grows like Ng(T). From [1], we know 
that N S (T) ~ C (7). We conclude that 

D p (Z l ) T ~ O (tf) . 

Finally, from (I14p . it is clear that 

which proves the theorem. 




3.2 The ^-variation estimates 

Similar to the quadratic variation estimate a 2 defined in ([7j), we define the p- variation 
estimates as the properly normalized total p- variation of the process: 

We will study the L2-error of this estimate in different scales. First, we need to define 
the constant C P (T). The natural choice would be to choose C P (T) so that E (a p ) = a p . 
So, 

C p (T) = hz((D p (Y^) T ) p ). 
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We need to compute E ((D p (Y 1 ' e ) T f) . From (Ej), we get that 

E ((A^ 1,£ )t) P ) = eVE ((p 9 {Z x ) t J) . 
Using (|T5|) . we get that 



E {{D p {Z') T y) = hmE [ £ \ W *» ~ W *\ P I > 

,t e e£ 5 ([o,T]) 



Note that for any p > 1, D V {Z X ) T < D X (Z X ) T , where E (/^(Z 1 )^) < +00. Thus, from 
the Dominated Convergence Theorem, the limit can come out of the expectation. To 
simplify our computations, from now on we will assume that Zq = Zf = 0. We have 
already observed that the random variables {(Wt /+1 — W te ), ti G Eg ([0, T]) , £ < Ng(T)~\ 
are independent and distributed like W TS where rg is the first time Z 2 crosses zero 
after t = 5, given that Zq = 0. Thus, 

E ( \ W ^-W te \ p \ =EN 5 (T)E\W TS \ p + E\W T -W tNsm \ p , 

\t t es 5 ([o,T\) J 

where Ng(T) is the number of zero-crossings of Z 2 in interval [0, T] that are distance 
5 apart from each other. First, we notice that 

E\W T5 \*> = E (E (\W Ts r\r 5 )) = -J=2ir f *±i) E (fo)*) . 

To compute E ((t$) p ), we note that rg can be written as rg = 5 + t(Z 2 ), where r{z) is 
the first zero-crossing of the process Z 2 given that it starts at z. For Z 2 an Ornstein- 
Uhlenbeck process, the p.d.f. of r{z) has been computed explicitly (see [T7]) and is 
given by 

tu , 2 \z\e~* ( z 2 e~ 2t " 

f(t,z) = —=- ^-rexp 



0F(l- e -a)§ V 1-e- 2 *, 

Since Zq = by assumption, Z| is a Gaussian random variable with zero mean and 
variance |(1 — e~ 2S ). Let us denote its p.d.f. by gg(z). It follows that the p.d.f. of 
t(Z 2 ) is given by 

, , , f°° , 4e*csch(5 + t) sinh(<5) sinh(t) 

where csch(t) = sin ^ t ^ and sinh(t) the hyperbolic sine. We write 

Kmi=E((7i)*) = hm-L !°° {5 + tfh 5 {t)dt = 

r (5 + t)Ph s (t) TTt , , 
= lim / ^ — KJ tH(t)dt. 
J y/StH(t) 
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where 



H(t) = { ] and J tH(t)dt = y/2. 



The function is increasing to £ 1+p as 5 4- an d thus by the dominated 

theorem we find that 

/•oo 

K p := / t p H(t)dt (18) 

3 

Notice that for t — > 0, behaves like and thus the integral K p is finite if and 
only if p > |. Also, for p = 1, 2 we find that = ^2 and K 2 = 2^/2 log 2. 

Now, we need to compute the limit of y/5¥,Ns(T) as 5 — > 0. We can use the results 
in [1] to get an upper and lower bound and show that Ng(T) behaves like O (^-^ 

However, we need to know the exact value of the limit. We proceed as follows: we 
write 

oo oo / n \ 

EN S (T) = ^P(iV,(T) >n) = J>(5>/<T (19) 

n=l n=l \ i=l / 

where rf — U — for tj e £g([0,T]) and i < Ng(T). Using (fTTj) we find that the 
Laplace transform of the distribution of Tg is 

ft (X) ~ e-»h (X) - jgj^lS j / A+ 1 A + 2 \ 
ff,(A)-e ^(A)- -j= _ == > e J, (20) 

where F\ (a, b, c, x) is the regularized hypergeometric function given by 

F 1 (a,6,c,x) = -i T f;^M^TT' and Wn = l[(d+k). 
mfo ^ k U to 

We find that for small d > 0, this behaves like 

o /ore A+l N i 

^(A) = l-^^|iVd+0(tf). (21) 

Since the rf's are i.i.d., the Laplace transform of the sum Xir=i ^ w ^ be ^ arL( ^ 
thus we write 



where £ 1 denotes the operator of the inverse Laplace transform. Substituting this 
back to (JT9~]1 . we get 

oo „y 



EiV,(T) = V / C-^Hg (X) n ](dt) 

n=l ^ 

oo 



n=l 

Vf ^1 ( 22 ) 
o I — tig [X) 
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Taking the limit inside the operator, we finally see that 
iV6E(N x (T)) = ^= T 

S— >c 



Finally, we note that since Z|. = by assumption, 

(T - t Ns{T) ) < 5 =► limE|Wr - W^ (T) | p = 0. 



For every p > 1, we set 



1 /p + 1\ % 



a„:=^2tr and c p :=-^if P . (24) 

Putting everything together, we find that 

¥,{{D p {Z l ) T ) P ) :=c p T (25) 

and consequently 

E ((D P (Y^) T ) P ) = e?a p c p ^ = e^ 2 a p c p T. 

Thus we set 

C P (T) := e p - 2 c p T. (26) 

By construction, the p-variation estimates a v defined in (fTB"]) are consistent, i.e. 
E (a p ) = a v . We now compute its square L2-error: 

El <m )-° = 

^ E (( D p( Yl ' e )T) 2P ) ~ (? 2P 

P 

f 2p /T 2p f / \ 2p 

E ((a^V 



C P (T) 2 



a 2 * 



<? P [4^({D P {Z')A 2P \-1) (27) 



To proceed, we need to compute the second moment of (D P (Z )t) ■ As with the 
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computation of the first moment, we write: 



E ((D^t)*) = limE Yl \ W ^~ W < 



t e \ 

t t ee s ([o,T\) 



\p 



limE \W T s\v+\W T -W tNs( , 



2 

\P 



71=1 



limE V \W T 



2 

IP 



5->0 I 



n=l 



where the last line comes from the fact that (T — ijv 5 (T)) < <5- To compute the above 
expectation, we write 



'N S (T) \~ / N S (T) 



MEW = E £ i^h^, 

n=l / \ m,n=l 



P 



= EiV 5 (T)E|J4y 2p + E (iV 5 (T) 2 - JVi(T)) (E|W T J P ) 2 
= EiV 5 (T)E|iy T ,| 2p + EiV 5 (r) 2 (E|WVj p ) 2 + C (v'd) 

where the last line follows from the fact that Ng(T) ~ O (^j^J an d ^|Wt«| p ~ 

o (ytj . it remains to compute the limit of 5KNg(T) 2 . Following a similar approach 
to the one before, we write 

oo oo / n 

EiV^T) 2 = J^(2n - 1)P (JV{(r) > n) = £(2n - 1)P ( ^ rf < T 

n=l n=l \i=l 



and 



in, \ oo /*T 

n =l \i=l / n=l ^° 

oo 

£^[^n^(A) n ](rft) 



n=l 

T 



y d 



Hg (A) 
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Taking the limit as 5 — > 0, we get 



Um5EN s (T) 2 = iim2 / C~H 5H * (A) J (dt) 



8^0 8^0 



H 5 (X) 




= ^ + (21og2)T 
Putting everything together, we get 

E {^{p p {Z 1 ) T ) 2p ^j = lim ^VSENs(T)E^^^- + 5EN S (T) 2 ^E-^^ 

= ^=a 2p K p + (2 log 2) 7^ (a p K f ) 2 

= Tc 2p + (T 2 + (4 log 2) T) (c p f 
= T 2 (c p ) 2 + T(c 2p +(41og2)(c p ) 2 ), 
where a p and c p are defined in (1241) and K p is defined in (I18|) . Finally, we get 

E^-o*) 2 = ^(^(^(c P ) 2 + ^(c 2p + (41og2)(c p ) 2 )) -l) 

= ^(^ + 41og(2))=a 2 ^( P ), (28) 

where £/(p) = + 4 log (2). This is an increasing function for p G [1,2] and 

41og2 =: E(l) < E(p) < E{2) := 10 log 2, Vp G [1,2]. 
We summarize our conclusions in the following 

Theorem 3.3. The L2-error of the estimator a p defined in (To)) is described by (H2J). 
At scale O (e a ), the error is of order O (e~ 



2-a 
2 



We see that the performance of the estimators a p is the same for all p > 1 and 
they outperform the a 2 estimator defined in (J7|). In terms of the constant E(p), the 
smaller the p, the smaller the error. However, there is a problem: except for scale 
(9(1) (a = 0), the normalizing constant C p depends on e, which will in general be 
unknown. We go on to define a new estimator that does not assume knowledge of e. 

3.3 Estimating the scale separation variable e 

Suppose that T < 1 and T = e a for some a > 0. We define the new estimator a p 
similar to a p , only use c p rather than C p as our normalization constant. Thus, we 
define 

a p = — (D p (Y 1,e ) T y , (29) 
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where c p is defined in ( 121T) . Then 



c 2 Prr 2 V / T 2 T \ f p rr p f T\ 

= — (7- (<*) + ? K + ( 41 °g 2 ) fe) )J - 2ct "— (<* ? J + " 2 " 

^(^ + ^(^ + (4.o g 2))- 2 (^) +1 ) (30) 

and by substituting T by e Q this becomes 

E (a p - a p f = o 2p ( e 2p+2a ~ A + e 2p+a ~ 2 E(p) - 2e p+a ~ 2 + l) . (31) 
Thus, we get the following behavior: 

(i) For p > 2 — a, the error is of order O (1). 

(ii) For p = 2 — a, the error is well-behaved and of order O {^~^~ 

(iii) For p < 2 — a and a < 2, the error explodes like O ( e 2 P+ 2a ~ 4 ). 

We conclude that the optimal estimator is <r 2 , since it does not assume knowledge 
of e and the estimators a p do not outperform it even for p = 2 — a (except that the 
constant E(p) is smaller). However, the estimators a p can be used to estimate the 
scale separation variable e. We set 



& 



p := arg mm 

Kp<2 



{a p y P -{a 2 Y 



and 

a := 2 — p. 

Then, we estimate e by 

e := T«. 
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